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Abstract 

Evolution of antibiotic resistance in microbes is frequently achieved by acquisition of spontaneous mutations during 
antimicrobial therapy. Here, we demonstrate that inactivation of a central transcriptional regulator of iron homeostasis 
(Fur) facilitates laboratory evolution of ciprofloxacin resistance in Escherichia coli. To decipher the underlying molecular 
mechanisms, we first performed a global transcriptome analysis and demonstrated that the set of genes regulated by Fur 
changes substantially in response to antibiotic treatment. We hypothesized that the impact of Fur on evolvability under 
antibiotic pressure is due to the elevated intracellular concentration of free iron and the consequent enhancement of 
oxidative damage-induced mutagenesis. In agreement with expectations, overexpression of iron storage proteins, inhi- 
bition of iron transport, or anaerobic conditions drastically suppressed the evolution of resistance, whereas inhibition of 
the SOS response-mediated mutagenesis had only a minor effect. Finally, we provide evidence that a cell permeable iron 
chelator inhibits the evolution of resistance. In sum, our work revealed the central role of iron metabolism in the de novo 
evolution of antibiotic resistance, a pattern that could influence the development of novel antimicrobial strategies. 

Key words: antibiotic resistance, bacterial evolution, SOS response. 



Introduction 

Understanding the molecular mechanisms promoting the 
evolution of antibiotic resistance is of crucial medical impor- 
tance. It has been suggested that prolonged exposure to sev- 
eral antibiotics induces changes in metabolism that stimulate 
intracellular accumulation of reactive oxygen species (ROS) 
(Dwyer et al. 2007; Kohanski et al. 2007, 2008). The impact 
of these reactive molecules on the evolution of antibiotic 
resistance is expected to be complex. On one hand, ROS 
formation has been implicated as a common step in antibi- 
otic-mediated cell death. As ROS interferes with normal phys- 
iological functions in respiring cells, these molecules are highly 
deleterious and damage several cellular subsystems (Imlay 
2003). ROS also leads to the accumulation of mutations 
either by directly damaging DNA and the nucleotide pool or 
through the activation of the SOS-mediated mutagenic DNA 
repair pathway (Cirz et al. 2005; Kohanski et al. 2010; Foti et al. 
2012). Thus, inducing the formation of ROS is potentially a 
double-edged sword, as it could be a major cause of both an- 
tibiotic-induced cell death and mutagenesis. The role of ROS in 
antibiotic killing is currently under intense debate (Keren et al. 
2013; Liu and Imlay 2013; Dwyer et al. 2014; Molina-Santiago 
and Ramos 2014) and is not a subject of this article. 

We focused instead on antibiotic-induced mutagenesis 
and examined the role of intracellular iron in this process. 
We investigated the hypothesis that perturbation of intracel- 
lular iron homeostasis has a key role in antibiotic-induced 



mutagenesis through the stress-induced formation of ROS. 
Indeed, numerous iron uptake-related genes are upregulated 
by bactericidal antibiotics, many of which are directly regu- 
lated by the ferric uptake regulator protein (Fur) (Dwyer et al. 
2007). Fur is a central transcriptional regulator of iron homeo- 
stasis and is a well-conserved protein present in many bacte- 
rial species. Upon binding Fe (II), Fur controls the expression 
of approximately 100 genes, many of which are involved in 
iron uptake and metabolism (Stojiljkovic et al. 1994; McHugh 
et al. 2003). Fur also has a major role in the regulation of acid 
shock response, oxidative stress response, toxins, and viru- 
lence factors (Escolar et al. 1999; Troxell and Hassan 2013). 
It is also considered a potential target for novel antimicrobial 
agents in Gram(— ) pathogenic bacteria (Horsburgh et al. 
2001; Mey et al. 2005). For all these reasons, we studied the 
impact of Fur on the evolution of antibiotic resistance. 

Although antibiotic resistance can emerge at very low an- 
tibiotic concentrations (Lopez and Blazquez 2009; Kohanski 
et al. 2010; Blazquez et al. 2012), studying development of 
resistance under lethal antibiotic stress is more relevant 
from a clinical point of view (Jumbe et al. 2003). Moreover, 
the evolution of resistance in this regime is expected to be far 
more complex, as antibiotic treatment may induce both 
mutagenesis and bacterial cell death. Therefore, our focus 
was on studying the evolution of resistance under lethal 
antibiotic pressure. We demonstrate that iron uptake, storage, 
and metabolism have a key role in stress-induced mutagenesis 
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during antibiotic stress in Escherichia coli K-12. Unexpectedly, 
deletion of a central regulator of iron homeostasis had no 
substantial effect on survival under toxic ciprofloxacin expo- 
sure, but facilitated the evolution of antibiotic resistance. Iron- 
overload-mediated mutagenesis during antibiotic stress was 
largely independent of the SOS response, indicating the pres- 
ence of parallel pathways. These results may indicate that 
blocking central components of the SOS pathway is insuffi- 
cient for eliminating the rise of resistant bacteria. Based on 
these results, we argue that oxidative damage in antibiotic- 
treated cells is a major driver of mutagenesis, especially when 
coupled with intracellular iron overload. Finally, we briefly 
discuss potential therapeutic consequences of our work: By 
binding unincorporated free iron, iron chelators can slow 
down the development of resistance. 

Results 

Inactivation of an Iron Homeostasis Regulator 
Promotes the Evolution of Antibiotic Resistance 
We investigated the evolution of resistance during a lethal 
ciprofloxacin exposure (6.25 times the minimal inhibitory 
concentration of wild-type [WT] £ coli). Ciprofloxacin is a 
widely used fluoroquinolone in the clinic, and its mechanism 
of action has been well studied (Barry et al. 1984; Eliopoulos 
et al. 1984). The dosage employed in this work is close to the 
mutant prevention concentration for £. coli K-12 that repre- 
sents a threshold above which the emergence of resistant 
mutants is not expected to occur (Olofsson et al. 2006; 
Olofsson and Cars 2007). All experiments were conducted 
in deep-well microtiter plates (at least 96 replicate 



populations per genotype). Plates were incubated for 5 
days, and the number of E. coli BW25113 resistant popula- 
tions was counted for each starting genotype. Ciprofloxacin is 
very toxic at such a high antibiotic concentration; therefore, 
the majority of the WT cultures with no resistant mutations 
became extinct by the end of the 5-day treatment (i.e., no 
growth observed after plating aliquots onto a nonselective 
medium). Evolution of resistance was observed only in 6-7% 
of the parallel populations established by WT cells (fig. 1A), 
most likely because adaptation to such a high dose demands 
rare or multiple specific mutations (Marcusson et al. 2009). 

To confirm the results of a previous genome-wide screen 
(Mehi et al. 2013) for the determinants of resistance evolu- 
tion, we first investigated the impact of the Fur regulator 
(Ferric uptake regulation) on the evolution of resistance in 
more detail. Afur populations showed an enhanced capacity 
to develop resistance: Approximately 50% of the indepen- 
dently evolving populations established by Afur cells were 
capable of acquiring resistance to ciprofloxacin, which is ap- 
proximately eight times more frequent than that observed in 
the isogenic WT strain (fig. 1A). This pattern is not due to an 
altered susceptibility of Afur to ciprofloxacin. Using an iden- 
tical experimental setup, viable cell numbers were determined 
by counting CFU (colony forming unit) values on Luria Broth 
(LB) agar plates. As expected, the bulk of the populations was 
killed in the first few hours of the treatment. Afur showed no 
significant increase in survival compared with the WT. If any- 
thing, a weak opposite trend was observed (fig. 1B); therefore, 
changes in ciprofloxacin sensitivity in Afur cannot explain its 
enhanced capacity to develop resistance. Daily monitoring of 
resistant clones in Afur populations during antibiotic 
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Fig. 1. Evolution of resistance against ciprofloxacin. Experiments were performed in the presence of 100ng/ml ciprofloxacin (6.25 times MIC of the 
WT). (A) Frequency of ciprofloxacin-resistant populations of the Afur genotype in comparison with the WT strain. ***P < 0.001 (chi-square test using 
the actual number of observations). Fractions are based on a total of n = 576 populations. Error bars represent confidence intervals of proportions 
(Wilson procedure). (B) Survival of Afur in the presence of ciprofloxacin compared with the WT strain. Error bars represent 95% confidence intervals. 
The mean of values is shown for n =3. n represents the number of repeated, independent experiments. 
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treatment showed that resistant mutants emerged de novo in 
the presence of antibiotic stress. Resistant clones were de- 
tected after only 2-3 days of incubation but showed fast 
growth after isolation and reinoculation into antibiotic-con- 
taining medium (data not shown). Based on these results, we 
conclude that deletion of Fur has a specific effect on antibiotic 
stress-induced mutagenesis. 

The Fur Regulon Is Plastic across Environments 
We examined the transcriptional responses of WT and A/ur 
cells with the aim of identifying differentially expressed genes 
that might influence the rate of resistance evolution. To 
achieve this goal, the transcription profiles were compared 
in the presence and absence of lethal ciprofloxacin stress 
(100ng/ml). Samples for microarray analysis were taken im- 
mediately before induction (time zero) and then at 60 min 
postinduction. Only genes that showed at least a 2-fold 
change in expression were considered further (see Materials 
and Methods). 

First, we compared the transcriptional profiles of the WT 
cells with those of A/ur grown in the absence of ciprofloxacin 
stress. A total of 125 genes were affected, of which 56 were 
induced and 69 repressed in A/ur (supplementary table SI, 
Supplementary Material online). Our data set was compared 
with results of a previous screen that aimed to identify Fur- 
regulated genes using microarrays (McHugh et al. 2003). 
In spite of the differences in experimental setup in the two 
works, there was a highly significant overlap in the sets of 
genes which changed their expression in A/ur (supplemen- 
tary fig. SI, Supplementary Material online). Gene ontology 
enrichment analysis (Boyle et al. 2004; Camon et al. 2004) 
revealed that genes involved in iron homeostasis, sidero- 
phore-mediated uptake, and enterobactin biosynthesis 
were upregulated in A/ur (supplementary table S2, Supple- 
mentary Material online). Remarkably, genes involved in an- 
aerobic respiration, cytochrome complex assembly, and 
members of the electron transport chain also had significantly 
changed expression levels. 

Ciprofloxacin caused a major reprogramming of gene ex- 
pression across the genome: 748 genes exhibited upregulation 
or downregulation (supplementary table SI, Supplementary 
Material online). In agreement with earlier works (Dwyer et al. 
2007), genes involved in DNA repair and SOS response were 
among the most highly upregulated genes (supplementary 
table S3, Supplementary Material online). Next, we compared 
the transcriptional responses of WT and A/ur cells in the 
presence/absence of ciprofloxacin. Interestingly, the Fur reg- 
ulon was substantially reorganized in response to ciprofloxa- 
cin stress (fig. 2). Altogether, 221 genes showed altered 
expression in A/ur under at least one of the two conditions 
(i.e., presence/absence of ciprofloxacin stress), of which only 
15.4% and 6.3% were up- and downregulated, respectively, 
in both conditions (fig. 2). In total, 96 genes were specifically 
affected in A/ur under ciprofloxacin stress, of which 57 were 
induced and 39 were repressed. To identify the most rele- 
vant cases, we applied a rigorous statistical procedure. The 
procedure identified genotype-environment gene expression 




Fig. 2. Low overlap in differently expressed genes upon fur deletion 
between ciprofloxacin-treated and untreated populations. Venn dia- 
grams showing the overlap in significantly upregulated (A) and down- 
regulated (B) ORFs upon fur deletion between ciprofloxacin-treated 
(CPR-treated) and untreated £ coli populations. The numbers indicated 
on the diagram refer to the number of genes with significantly altered 
expression level in one or both conditions. 



interactions at the level of individual genes (see Materials and 
Methods) and revealed several genes that show strong up- or 
downregulation in A/ur only under ciprofloxacin stress 
(table 1). For example, the mRNA level of SufA, a scaffold 
protein involved in Fe-S cluster assembly, remained largely 
unaltered in response to ciprofloxacin stress or deletion of fur. 
In contrast, the expression level of the same gene showed an 
8-fold upregulation in A/ur under ciprofloxacin stress com- 
pared with WT cells under no stress (table 1). Other gene 
products included a member of the major iron transport 
membrane system (ExbB) and two enzymes (EntA, EntC) 
involved in enterobactin biosynthesis. Based on these pat- 
terns, we hypothesized that the impact of Fur on antibiotic 
resistance is specifically linked to its effect on iron uptake. This 
hypothesis was investigated further through in-depth genetic 
analyses. 

Impact of Iron Uptake and Storage on Resistance 
Evolution 

As the Fur protein activates iron-storage proteins and 
represses iron uptake, inactivation of the corresponding 
gene yields enhanced intracellular concentration of reactive 
ferrous iron (Abdul-Tehran i et al. 1999). We first investigated 
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Table 1. 


Genes Showing Significantly Altered Transcriptional Response in Afur in 


the Presence of Ciprofloxacin Compared with What Is Expected 


Based on 


an Independent Effect of Afur and Ciprofloxacin. 








Gene 


Function of the Encoded Protein 




mRNA Level Relative to Untreated WT 






Afur 


CPR 


Afur + CPR 


exbB 


Subunit in the energy-transducing TonB-ExbB-ExbD complex 


1.43 


1.1 


6.64 


efeU 


Hypothetical protein, subunit of ferrous iron transporter 


1.05 


■ 


1 72 


mntH 


Mn 2+ /Fe 2+ : H+ symporter 


1.05 


U.o 1 


1 .DO 


yq/'H 


Siderophore interacting, NADPH-dependent ferric reductase 


1.62 


n 


z./o 


entC 


Isochorismate synthase involved in enterobactin biosynthesis 


3.39 


n £q 
u.oy 




sufA 


Fe-S cluster assembly scaffold protein 


1.7 


1.17 


8.25 


entA 


2,3-Dihydro-2,3-dihydroxybenzoate dehydrogenase involved in 


1.81 


0.96 


5.89 




enterobactin biosynthesis 








ybbB 


Selenophosphate-dependent tRNA 2-selenouridine synthase 


1.01 


1.71 


1.23 


ptrB 


Protease II 


1.1 


1.24 


1.1 


lolB 


OM lipoprotein required for localization of lipoproteins 


1.16 


2.53 


1.8 


IpxL 


Lauroyl acyltransferase 


1.21 


1.47 


1.22 


napF 


Predicted 3Fe-4S iron sulfur protein, stimulates the activity of 


0.91 


3.85 


1.73 




periplasmic nitrate reductase (Nap) 









Note. — Only those ORFs are listed which showed a significant interaction term in the two-way analysis of variance test after controlling the false discovery rate at 10%. Higher 
(lower) than expected transcript levels in A/ur in the presence of ciprofloxacin (CPR) are indicated in white (gray) background color. 



the impact of iron uptake on the evolution of antibiotic 
resistance. TonB is the main energy-transducing subunit of 
numerous, high-affinity outer membrane iron transporter 
complexes; therefore, ferric ion transport can partially be in- 
hibited by suppressing the activity of this gene (Andrews et al. 
2003). In contrast to a prior study in Pseudomonas aeruginosa 
(Yeom et al. 2010), inactivation of TonB had no significant 
effect on survival under toxic ciprofloxacin exposure in £ coli 
(fig. 3A). However, AfurAtonB double mutants showed a 
markedly reduced capacity to develop resistance compared 
with that of Afur (fig. 3B). Thus, accelerated evolution 
of resistance in Afur is contingent upon intracellular 
accumulation of free iron and it can be suppressed by 
inactivating TonB. Next, we tested the impact of intracellular 
iron sequestration on the development of resistance thr- 
ough overproduction of iron-storage proteins. FtnA and Bfr 
encode the distantly related iron-storage proteins ferritin 
and bacterioferritin, respectively. The plasmids carrying 
these genes were introduced individually into the Afur 
background. Overproduction of any of these iron-storage 
proteins suppressed resistance evolution (fig. 3C). Taken 
together, both extracellular iron uptake and internal iron 
storage have crucial roles in Fur-mediated evolution of 
resistance. 

Links between Iron Overload and Antibiotic 
Stress-Induced Mutagenesis 

It has been suggested that bactericidal antibiotics induce mu- 
tagenesis by stimulating the formation of ROS (superoxide, 
peroxide, and hydroxyl radicals) that results from a perturba- 
tion of the citric acid cycle (Dwyer et al. 2007; Kohanski et al. 
2007). The key assumption of the hypothesis is that bacteri- 
cidal antibiotics (such as ciprofloxacin) promote acceleration 
of the Fenton reaction: Fe 2+ + H 2 O z Fe 3+ HO* + H 2 0. HO* 



is a highly reactive oxidant that promotes mutations by di- 
rectly damaging the nucleotide pool and DNA (Farr et al. 
1986; Touati et al. 1995; Nunoshiba et al. 1999, 2002). As 
both iron uptake and storage influence the evolution of re- 
sistance, we suggest that the rate of the Fenton reaction 
during antibiotic exposure may be limited by the intracellular 
availability of ferrous iron. This hypothesis was investigated in 
detail. 

As it has remained a contentious issue (Liu and Imlay 
2013), we first measured the accumulation of ROS during 
antibiotic stress. Dihydrorhodamine 123 dye (DHR) is suitable 
for use as an indicator of oxidative stress (Henderson and 
Chappel 1993; Yeom et al. 2010). Compared with untreated 
controls, ciprofloxacin-treated populations oxidized the dye 
more efficiently (fig. 4A). However, as the dye is not comple- 
tely selective, the identity of the accumulated oxidants 
remained uncertain in this experiment. To directly test 
the contribution of superoxide and H 2 0 2 formation under 
antibiotic treatment, we used GFP reporter plasmids, 
and measured changes in promoter activities of genes in re- 
sponse to antibiotic stress (Zaslaver et al. 2006). We moni- 
tored promoter activities of 1) OxyR, a transcription factor 
induced by H 2 0 2 stress, 2) members of the OxyR regulon with 
specific roles in H 2 0 2 scavenging (KatG, AhpC), and 3) the 
SoxS dual regulator, which is generally activated by superoxide 
and redox cycling compounds (Farr and Kogoma 1991). 
The 2- to 5-fold increases in promoter activities were 
observed in all cases (fig. 4B-D). These results suggest that 
ciprofloxacin promotes a mild but significant accumulation 
of ROS. 

We next tested whether the enhanced capacity to de- 
velop resistance in Afur is linked to respiration and 
subsequent generation of superoxide molecules. The man- 
ganese (SodA) and iron (SodB) superoxide dismutases are 
responsible for detoxification of superoxide in £ coli, and 
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Fig. 3. Impact of iron uptake, iron storage, and inactivation of error prone polymerases on ciprofloxacin-resistant evolution. (A) Effect of tonB deletion 
on survival under ciprofloxacin exposure in WT and Afur genetic backgrounds. The mean of values is shown for n = 8. (B) Effect of impaired ferric iron 
uptake due to deletion of tonB gene on the fractions of ciprofloxacin-resistant populations in the case of Afur and AsodAB. Fractions are based on a 
total of n = 96 (wt vs. ctonB), n = 288 (Afur vs. Afur AtonB), and n = 192 (AsodAB vs. AsodABctonB) populations of each genotype. Error bars represent 
confidence intervals of proportions (Wilson procedure). (C) Effect of overexpressing mutS and genes involved in iron storage (bfr,ftnA) on the fractions 
of ciprofloxacin-resistant populations in the case of Afur. Fractions are based on a total of n = 96 populations of each genotype. Error bars represent 
confidence intervals of proportions (Wilson procedure). Genes of interests were overexpressed from a high-copy number plasmid (pZE31) in a Afur 
genetic background. Control, Afur harboring empty pZE31 plasmid; ++, overexpression of a given gene from the pZE31 plasmid. (D) Effect of the 
inactivation of error prone DNA polymerases PollV (dinB) and PolV (umuDC) on the fractions of ciprofloxacin-resistant populations in the case of Afur. 
Fractions are based on a total of n = 288 populations of each genotypes. Error bars represent confidence intervals of proportions (Wilson procedure). 
***P < 0.001; ns, P > 0.05 (chi-square test using the actual number of observations). 



both genes are part of the Fur regulon (Farr and Kogoma 
1991; McHugh et al. 2003). As expected, deletion of both 
genes (AsodAB) enhances mutation rate (Nunoshiba et al. 
1999, 2002) and the consequent evolution of ciprofloxacin 
resistance (fig. 3B). More remarkably, simultaneous inacti- 
vation of tonB suppressed the evolution of resistance in 
AsodAB (fig. 3B). Therefore, the mutagenic effect of en- 
hanced superoxide formation was contingent upon iron 
uptake and probably relies on Fenton chemistry, similarly 
to that observed in Afur. As a further support, we exam- 
ined the evolution of ciprofloxacin resistance under 



anaerobic growth conditions. The frequency of resistance 
was considerably reduced in parallel evolving Afur and 
AsodAB populations in the absence of oxygen (fig. 5A). 
This result is unlikely to reflect differences in antibiotic 
selection pressure between aerobic and anaerobic condi- 
tions, as ciprofloxacin had a substantial killing effect inde- 
pendent of oxygen availability (Malik et al. 2007; Wang 
et al. 2010) (fig. 5A). We conclude that the mechanisms 
underlying enhanced evolvability in Afur and AsodAB are 
linked and are most likely associated with antibiotic- 
induced oxidative stress. 
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Fig. 4. Direct measurement of oxidative stress induced by ciprofloxacin treatment. (A) Accumulation of ROS following exposure to 100 ng/ml 
ciprofloxacin in Afur, AsodAB, and WT genotypes. DHR fluorescence was measured 24 h after addition of ciprofloxacin. Final fluorescence intensities 
were normalized to cell densities (OD600) and are shown relative to the untreated WT control. The mean of values is shown for n = 3. (B-D) Induction 
of promoters related to oxidative stress: oxyR, OxyR-regulated katC, ahpC, and soxS following exposure to 100 ng/ml ciprofloxacin. GFP intensities, 
indicating promoter activities, were measured 24 h after addition of ciprofloxacin. Relative promoter activities are shown; linear GFP signal intensities 
were normalized with GFP intensities of the untreated WT (B, C) or with GFP intensities of the untreated wild type bearing a promoterless construct 
(control in panel [D]). Error bars represent 95% confidence intervals. *P< 0.05; **P< 0.01; ***P< 0.001; ns, P > 0.05 (two-sample t-test). The mean of 
values is shown for n = 2. n represents the number of repeated, independent experiments. For more details, see Materials and Methods. 



Finally, we tested the impact of iron chelators on the 
evolution of resistance. Iron chelators interfere with the 
Fenton reaction, and thereby they inhibit generation of 
hydroxyl radicals. New sets of experiments were initiated 
with the same founding strains and antibiotic settings, with 
the cell permeable iron chelator o-phenantroline also added 
to the medium. Importantly, cotreatment with subinhibitory 
concentration of o-phenantroline had no measurable ef- 
fect on ciprofloxacin minimum inhibitory concentration 
(data not shown). At the same time however, it dramatically 



reduced the frequency of resistant populations in both wild 
type and Afur populations (fig. SB). 

Iron-Overload-Mediated Mutagenesis Is Partly 
Independent of Error-Prone DNA Polymerases and 
Mismatch Repair 

To investigate the underlying molecular mechanisms of mu- 
tagenesis in more detail, we assessed the involvement of 
error-prone DNA polymerases PollV (dinB) and PolV 
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(umuDC). These enzymes have an important role in SOS-de- 
pendent stress-induced mutagenesis, and the DinB protein is 
directly involved in the incorporation of oxidized guanine into 
DNA (Friedberg et al. 2002; Foti et al. 2012). We constructed 
the AdinBAfur and AdinBAumuDCAfur mutants, and ini- 
tiated new sets of experiments with these strains. AfurAdinB 
retained the enhanced capacity to develop resistance, 
whereas simultaneous deletion of PollV and PolV caused 
minor, but significant reduction in the frequency of resistant 
populations (fig. 3D). Taken together, these results indicate 
that the enhanced mutagenesis in Afur is partly independent 
of the activity of the error-prone polymerases. Based on prior 
work, it has been suggested that the efficiency of the mis- 
match repair system (MMR) is decreased during stationary 
phase and certain antibiotic stresses (/3-lactam antibiotics) 
due to the transient depletion of the mismatch repair detect- 
ing and binding MutS protein (Gutierrez et al. 201 3). For these 
reasons, we tested whether transient depletion of MutS is 
responsible for the elevated rate of resistance evolution in 
Afur. This was performed by introducing a mutS-carrying 
plasmid into Afur. The elevated capacity of Afur to develop 
resistance was not abolished by overproducing MutS (fig. 3C). 
This result is in good agreement with the previously reported 
inability of the MMR to recognize and repair ROS induced, 
oxidative DNA lesions (Schaaper et al. 1989; Fowler et al. 2003; 
Wyrzykowski and Volkert 2003). Moreover, a recent work 
identified a significant decrease in mutation rates in MMR 
deficient, laboratory evolved £ coli lines associated with en- 
hanced ROS detoxification (Turrientes et al. 2013). The weak 
effect of deactivated error-prone DNA repair in AdinBAfur 
and AdinBAumuDCAfur mutants indicates that the major- 
ity of resistance mutations is formed by direct oxidation of 
DNA and not by incorporation of oxidized nucleotides 



through error-prone DNA polymerases (Foti et al. 2012). 
Indeed, we failed to find significant differences in the spectra 
of accumulated adaptive mutations between WT and Afur 
populations (supplementary fig. S2, Supplementary Material 
online). 

Discussion 

In times of extreme antibiotic stress, bacteria benefit from an 
enhanced mutation rate, as this can provide mutations that 
protect the population from death. Microbes have a broad 
range of molecular toolkits that promote adaptation, includ- 
ing stress-induced mutagenesis pathways (Bjedov et al. 2003; 
Kohanski et al. 2010; MacLean et al. 2013). Here, we demon- 
strate that exposure to a toxic dosage of ciprofloxacin 
boosts bacterial evolution partly through an iron-overload- 
dependent mutagenesis pathway (fig. 6). The main findings of 
our work are as follows: First, inactivation of a central regula- 
tor of iron homeostasis (Fur) promotes laboratory evolution 
of ciprofloxacin resistance. Second, intracellular iron uptake 
and storage plays a central role in this process. Third, super- 
oxide formation correlates with the capacity to develop 
resistance, and this capacity is dependent on intracellular 
iron availability. Fourth, many of the mutations delivering 
resistance in Afur are not the result of stress-induced muta- 
genesis through error-prone polymerases IV and V. One pos- 
sibility is that the main replicative polymerase (Pol III) is 
responsible for the high level of mutations through its re- 
cently identified capacity of incorporating oxidized nucleo- 
tides (Yamada et al. 2012). 

In summary, it appears that iron homeostasis influences 
the evolution of ciprofloxacin resistance. This enhanced evo- 
lutionary capacity is accompanied by only a minor reduction 
in population survival under antibiotic stress (fig. 1B). 
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Fig. 6. Schematic model of iron-overload-dependent oxidative mutagenesis in the Afur mutant, exposed to a toxic dose of ciprofloxacin. The primary 
antibiotic-target (ciprofloxacin-DNA gyrase) interaction induces physiological changes in the bacterial cell, including perturbations in metabolism and 
respiration. The electron transport chain becomes hyperactivated, which stimulates superoxide (0 2 ) formation. In the absence of Fur, the reduced 
activity of superoxide dismutases (SodA and SodB) further increases the intracellular level of superoxide. Superoxide damages Fe-S clusters of the 
proteins leading to the release of Fenton reactive ferrous iron (Fe 2+ ). The constitutive siderophore-mediated ferric iron (Fe 3+ ) uptake and decreased iron 
storage activity, together with the antibiotic-mediated iron release leads to an intracellular iron overload. Oxidation of the ferrous iron in the Fenton 
reaction results in the formation of the highly reactive hydroxyl radical (OH ), which may damage nucleotides and DNA, leading to the formation of 
mutations. Numbers represent intervening points by which antibiotic resistance development could be significantly reduced: 1) Inhibition of side- 
rophore-mediated iron uptake by inactivation of tonB, 2) enhancement of iron storage by overexpressing FtnAand Bfr iron storage proteins, 3) chelation 
of intracellular free iron using a cell permeable iron chelator o-phenantroline, and 4) applying anaerobic conditions. 



Why is this so? We speculate that oxidative damage in 
antibiotic-treated cells is a major driver of mutagenesis 
(especially when coupled with intracellular iron overload), 
but its contribution to antibiotic killing is relatively minor 
in our experimental settings. This explanation is in accor- 
dance with the existence of two, partly independent path- 
ways mediating ciprofloxacin-induced cell death, where only 
one of them is driven by hydroxyl radicals (Wang et al. 2010). 
We emphasize that these results were reached at an antibiotic 
dosage only six times above the WT MIC. It remains to be 
elucidated whether killing significantly prevails over mutagen- 
esis at very high antibiotic concentrations. Future work 
should also explore the prevalence of iron metabolism affect- 
ing mutations in drug-resistant pathogens. 

Finally, our work has potential implications for the devel- 
opment of future antimicrobial strategies. By binding unin- 
corporated free iron, iron chelators can slow down the 
development of resistance (fig. 5B). More generally, targeting 
bacterial iron uptake (Ballouche et al. 2009) and iron metab- 
olism using novel compounds is a double-edged sword, 
as iron appears to be critical for both bacterial growth and 
mutagenesis. 

Materials and Methods 

Bacterial Strains and Growth Conditions 

Unless otherwise indicated, this work employed £ coli 

BW25113 (K-12) and corresponding deletion mutants 



selected from the KEIO collection (see table 2 for strain list; 
Baba et al. 2006). Multiple, markerless deletions were achieved 
using subsequent P1 transductions (Green and Sambrook 
2012) and removal of selection markers. For this purpose, 
plasmid borne (pFT-A) FLP recombinase was used in order 
to excise kanamycin resistance cassettes (Posfai et al. 1997). 
All strains were retested by polymerase chain reaction (PCR) 
for the presence of multiple deletions and the loss of kana- 
mycin cassette. GFP reporter plasmids were purified from a 
comprehensive promoter library (Zaslaver et al. 2006). All 
plasmids were transferred into WT and different mutants 
of interest (table 2) by electroporation. In all experiments 
involving expression vectors and reporter constructs, 
we used the empty plasmids as negative controls. Overexpres- 
sion of Bfr, FtnA, and MutS proteins was achieved by subclon- 
ing of PCR-amplified ORF sequences from the corresponding 
ASKA library plasmids (Kitagawa et al. 2006) into the pZE31 
vector (Lutz and Bujard 1997). HIS-tags that were originally 
introduced into coding sequences were removed by recon- 
structing native, ATG start-codons using specific cloning 
oligonucleotides (supplementary table S4, Supplementary 
Material online). In the case of all primer-pairs, a Hindlll, 5'- 
overhanging recognition site was utilized for cloning pur- 
poses. The resulting overexpression vectors were introduced 
into Afur and WT genetic background lacking the cognate 
TetR repressor protein (Lutz and Bujard 1997) allowing a fully 
derepressed state of expression in these strains. Unless 
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Table 2. The List of Strains and Plasmids Used in This Work. 



Genotype or Relevant Characteristics 



Source 



Strains 
BW25113 



AtonB 
Afur 

A/ur AtonB 
AfurAdinB 
AfurAdinBAumuDC 
AsodAB 
AsodABAtonB 
Plasmids 
pZE31 
pZE31 Bfr 
pZE31 FtnA 
pZE31 MutS 
pUA66 
pUA66 oxyR 
pUA66 katG 
pUA66 _ahpC 
pUA66 soxS 



rrnB3 AlacZ4787 hsdR514 A(araBAD)567 

A(rhaBAD)568 rph-1 
Same as 

BW25113, AtonB, KanR 

BW25113, Afur (removed Kan cassette) 

Afur, AtonB 

Afur, AdinB 

AfurAdinB, AumuD, AumuC 
BW25113, AsodA, AsodB 
AsodAB, AtonB 

High copy, ColEI, PLtetO-1, CmR 

pZE31, inserted ORF of Bfr 

pZE31, inserted ORF of FtnA 

pZE31, inserted ORF of MutS 

Low copy, pSOOl Ori, KanR 

pUA66, inserted promoter region of oxyR 

pUA66, inserted promoter region of katG 

pUA66, inserted promoter region of ahpC 

pUA66, inserted promoter region of soxS 



Baba et al. (2006) 



Baba et al. (2006); this work 
This work 
This work 
This work 

Baba et al. (2006); this work 
This work 

Lutz and Bujard (1997) 
This work 
This work 
This work 

Zaslaver et al. (2006) 
Zaslaver et al. (2006) 
Zaslaver et al. (2006) 
Zaslaver et al. (2006) 
Zaslaver et al. (2006) 



Note. — KanR, kanamycin resistance; CmR, chloramphenicol resistance. 



otherwise indicated, all experiments were carried out in LB or 
filter-sterilized LB (FACS experiments) at 37 °C (Green and 
Sambrook 2012). For the experiments involving the iron che- 
lator, growth medium was supplemented with o-phenan- 
throline (Sigma) at subinhibitory dosage (50 uM). Anaerobic 
experiments were carried out in a Bactron X anaerobic work- 
station under anaerobic atmosphere (95% N 2 /5% H 2 with 
palladium catalyst). Final concentrations of kanamycin (re- 
porter constructs) and chloramphenicol (overexpression con- 
structs) were 50 and 20(ig/ml, respectively. Antibiotics and 
chemicals were purchased from Sigma and Becton Dickinson 
(Bacto). 

Measuring Intrinsic Antibiotic Susceptibility and 
Evolution of Resistance 

Minimum inhibitory concentrations (MIC) for the WT and 
different mutant strains were measured by a standard mi- 
crodilution method in 96-well microtiter plates. We estab- 
lished a short term, high-throughput screen to study de novo 
evolution of resistance against ciprofloxacin (6.25x MIC; 
Mehi et al. 2013). Overnight cultures of each genotype were 
grown in deep-well plates. From each of these cultures, large 
populations (~10 8 cells) were transferred into 350 uJ fresh LB 
medium supplemented with a single antibiotic. In total, there 
were at least 96 replicate evolving populations per genotype 
and treatment. Deep-well plates were covered with sandwich 
covers (Enzyscreen) and incubated in a gyratory shaker 
(37 °C, 320 rpm). Following 5 days of incubation, 1—2 p.1 ali- 
quots from each well were transferred to LB agar plates with 
or without supplemented antibiotic. Growing patches of 



surviving populations (LB) and resistant populations 
(LB + antibiotic) were counted for each genotype and treat- 
ment. Due to slow growth under anaerobic conditions, cells 
were incubated for 3 days without shaking prior to transfer 
into antibiotic-supplemented medium, and adaptation 
against ciprofloxacin lasted for 11 days. Following anaerobic 
adaptation, screening for resistant populations was carried 
out as described above. 

Cell Viability Assays 

Cell survival was measured by determining changes in popu- 
lation size upon exposure to ciprofloxacin (100 ng/ml) in WT 
and various mutant genetic backgrounds. The experiments 
were conducted in 96 deep-well masterblock plates (350 ul LB 
medium supplemented with the antibiotic). Each well initially 
contained approximately 10 8 cells. During antibiotic treat- 
ment, samples were taken at multiple time points (0, 1, 3, 
and 6h postexposure) from 4-4 parallel populations per 
strain. The total number of viable cells (CFU) was determined 
by plating and counting diluted cultures on LB agar plates. By 
supplementing agar plates with ciprofloxacin (100 ng/ml), we 
could confirm that the number of resistant cells remained 
undetectable during the entire course of these experiments. 
Additionally, we measured survival under anaerobic condi- 
tions using the same technique as described above. 

Promoter Activity Measurements 
To measure changes in promoter activity in response to cip- 
rofloxacin stress, we used selected constructs from a compre- 
hensive GFP transcription reporter library and transferred 
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them into different genetic backgrounds (table 2) (Zaslaver 
et al. 2006). Each reporter strain bears a low-copy plasmid 
with a promoter of interest controlling the expression of a fast 
folding GFP and carries a kanamycin resistance cassette. The 
list of investigated genes includes soxS, oxyR, ahpC, and katG. 
Overnight cultures of reporter strains were preincubated in 
microtiter plates for 24 h in 100 ul filter-sterilized LB supple- 
mented with 50ug/ml kanamycin. Each strain was repre- 
sented by eight parallel populations per plate. Fifty 
microliters of each culture was transferred to fresh LB. In 
the fresh medium, the final concentration of ciprofloxacin 
was set to lOOng/ml. In total, there were four replicate pop- 
ulations per genotype and treatment (four untreated and 
four treated). Flow-cytometric measurements were per- 
formed using a Millipore Guava 8HT instrument following 
24 h of postexposure ciprofloxacin treatment. Mean, linear 
GFP signal intensities were calculated from 15,000 events 
per population and were normalized to untreated, WT 
signal intensities. Events lacking fluorescence signal in green 
or red (autofluorescence) were gated out from the analysis. 

Direct Measurement of ROS Production 
Evaluation of ROS production during ciprofloxacin treatment 
was conducted by a set of fluorescence measurements in a 
BioTek Synergy2 plate-reader. As a redox sensitive fluorescent 
dye, we applied dyhidrorhodamine123 (Sigma) in the pres- 
ence and absence of 100ng/ml ciprofloxacin (Yeom et al. 
2010). Initially, mutant (Afur, AsodAB) and WT strains 
were inoculated from single colonies and incubated for 24 h 
in Erlenmeyer flasks at 37 °C Antibiotic treatment and fluo- 
rescence measurements were carried out in microtiter plates 
covered with sandwich covers (Enzyscreen) after adjusting 
initial cell densities to OD 600 = 1.5 in each population. 
Control and treated populations (9-9 replicates per strain) 
were preincubated for 3 h prior to the addition of the fluo- 
rescent dye (10).ig/ml) in a gyratory shaker at 300 rpm and 
37 °C Following additional 24 h of incubation, fluorescence 
induction was measured in a Synergy2 fluorescent plate- 
reader (with the combination of 485/20 nm-excitation and 
528/20 nm-emission band-pass filters). In all cases, final fluo- 
rescence intensities were normalized to cell density (OD600). 
We calculated and corrected within-plate effects by measur- 
ing fluorescence intensities of untreated WT populations 
under the same conditions. 

Microarray Analysis 

We compared the microarray-determined mRNA profiles 
(Affymetrix E. coli Genome 2.0 Array) of Afur and WT geno- 
types in response to ciprofloxacin treatment (100 ng/ml) with 
those of untreated cultures. We used two biological replicates 
per genotype. Cultures were started from single colonies and 
grown overnight in Erlenmeyer flasks (25 ml LB). Samples for 
microarray analyses were taken before ciprofloxacin treat- 
ment (tO) and then at 1 h posttreatment (t1). Cell concen- 
tration was adjusted to 10 9 cells/ml and cultures were treated 
with ciprofloxacin in Erlenmeyer flasks using 20 ml fresh LB. 
Approximately 10 9 cells were taken at each time point as 



samples for RNA isolation. In order to stabilize RNA, RNA 
Bacteria Protect Reagent (Qiagen) was added to the samples 
according to the instructions of the manufacturer and 
then samples were stored overnight at —80 °C Total RNA 
was isolated by using RNeasy Mini Kit (Qiagen) according to 
the manufacturer's instructions. Sample concentration was 
estimated using the NanoDrop 1000 (Thermoscientific) 
spectrophotometer. 

We processed and analyzed raw gene expression data-files, 
using the affy package 1.38.1 (Gautier et al. 2004) in the R 3.0.2 
(R Core Team 2013) programming environment. Probesets 
were annotated based on the ecoli2.db 2.9.0 package 
(Gentleman et al. 2004). In order to reduce chip-to-chip var- 
iation, first we RNA normalized the intensity values of the 
probesets with assigned E. coli locus IDs and subsequently 
calculated expression fold changes. Robust Multi-array 
Average (RMA)-normalized data for each gene and each con- 
dition can be found in supplementary table S5, 
Supplementary Material online. Also full, RMA normalized 
data can be accessed at NCBI GEO (GSE55662). We chose 
the threshold value for up- and downregulated genes as 2- 
fold, which yielded a highly significant overlap with an earlier 
gene expression study of the fur regulon (McHugh et al. 2003) 
(hypergeometric test, P = 3.5x 10~ 25 and P=5x 10~ 4 for 
up- and downregulated genes, respectively; see also supple- 
mentary fig. SI, Supplementary Material online). In order to 
test the interaction between ciprofloxacin treatment and 
gene knockout, we employed an analysis of variance test 
with 10% false discovery rate correction. 

Supplementary Material 

Supplementary tables S1-S5 and figures S1 and S2 are avail- 
able at Molecular Biology and Evolution online (http://www. 
mbe.oxfordjournals.org/). 
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